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Abstract 

We find that Gamma-ray bursts (GRBs) at the 3B catalog are corre- 
lated (at 95% confidence level) with Abell clusters. This is the first known 
correlation of GRBs with any other astronomical population. It confirms 
the cosmological origin of GRBs. Comparison of the rich clusters auto- 
correlation with the cross-correlation found here suggests that ~ 26 ± 15 % 
of an accurate (S < 2.3°) position GRBs sub-sample members are located 
within 600 h -1 Mpc. 

1 Introduction 

One of the major obstacles in GRBs research is the lack of any association of GRBs 
with any other astronomical population. The recent observations of the BATSE experiment 
on the COMPTON-GRO observatory suggest that GRBs are cosmological (Meegan et al. 
1992.) The GRBs distribution appears to be isotropic and there is a paucity of weak bursts. 
Both facts are naturally explained by a cosmological distribution. The observed peak flux 
distribution in the BATSE catalog agrees well with a theoretical cosmological peak flux 
distribution (Piran 1992; Mao & Paczyhski 1992; Dermer 1992; Wickramasinghe et al. 1993; 
Cohen & Piran 1995.) These facts can, however, be accommodated by some extended 
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Galactic halo models and there is an ongoing debate on the galactic or extra-galactic origin 
of GRBs (Paczyhski 1995; Lamb 1995.) 

We have discovered a positive cross-correlation signal between the GRBs distribution 
in the 3B catalog and Abell clusters. This confirms our previous result that indicated the 
existence of such a correlation in the BATSE 2B catalog (Cohen, Kolatt, & Piran 1994.) 
This correlation demonstrates the cosmological origin of GRBs. It also enables us to estimate 
directly the distance scale to GRBs. We did not find any significant auto-correlation of GRBs 
in the 3B catalog. Since GRBs are distributed over cosmological distances (z « 1) one should 
not expect a cross-correlation with galaxy catalogs that consist of relatively nearby objects. 
Indeed we did not find any cross-correlation with IRAS galaxies neither with galaxies in any 
other optical catalog. However, The volume limited sample of Abell clusters (Abell 1958; 
Abell, Corwin, & Olwin 1989) extends to larger distances and may lead to a detectable 
cross-correlation. 

2 The Data 

The data consist of 3616 Abell clusters of richness class R > in Galactic latitudes 
\b\ > 30°. The cutoff is meant to reduce statistical noise emerging from low values of the 
cluster selection function cf> c i (Scaramella et al. 1991.) The clusters show a strong auto- 
correlation signal (Bahcall & Soneira 1983; Batuski et al. 1989; Peacock & West 1992), and 
their typical redshift values 0.15 — 0.2 . There are 549 GRBs in the 3B catalog 

with \b\ > 30°. The GRBs selection function, (f> gr b } is anisotropic due to different exposure 
durations in different directions. It can be derived from the exposure duration ratio as 
function of declination and has up to 20% variation across the sky. The position of the 



bursts is uncertain, with an angular error for the j th burst given by 8j = J 8 2 + 8 2 si 



where 6 sys = 1.6° is a systematic error and S sta t,j is a statistical error that varies from one 
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burst to another. The statistical error, 8 sta t,j^ decreases with the burst strength. For a weak 
burst 8 sta t,j can be quite large (up to 20°) and it might be meaningless to include objects 
with such poorly known positions in our analysis. Hence we have constructed an "accurate" 
sub-class of the GRBs sample, that contains 136 GRBs with a good positional accuracy, 
namely for which the positional error 8 3 satisfies: 8 3 < 8 max = \f28 sys . Clearly, because of 
the nature of this sub-class it contains stronger bursts on average. 



3 Cross Correlation Estimate 



The cross-correlation of the Abell clusters and the anisotropic sky coverage of the 
BATSE data make it difficult to estimate the cross-correlation between these two data sets 
in the usual way. We do like however to start up with the conventional cross-correlation 
evaluation in order to obtain its model independent values. Let N gr b_ c i(6) be the number of 

pairs of GRBs and clusters separated by the angle where for the i th bin < < 0i and 
let N gr b_ po (6) be the number of pairs of GRBs and randomly distributed particles subject to 
the cluster selection function </> c /, then the angular cross-correlation function, w(6), is defined 



as 



l + w(9)= ^ rb - d 2 —> (1) 
N grh - V0 [V) n d 

where n c \ and n po are the number densities of the clusters and Poisson particles respectively. 
In order to obtain a bias free estimate of the correlation, we corrected for the GRBs selection 
function by assigning a weight to each pair count inversely proportional to <f> grb . 

The main advantage of the cross-correlation assessment is that it allows to reduce the noise 
level appreciably in comparison to the GRBs auto-correlation. In the limit where n po ^> n c i } 
the Poisson error in w(0) is (Aw) p = (1 + w)(N gr b- c i)~ 1 ^' 2 for each bin. The overall error 
in the correlation value is larger than the Poisson error due to the cluster auto-correlation, 
the possible GRBs auto-correlation, and the assigned weights. Considering the cluster auto- 
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Figure 1: The probability distribution for random artificial GRBs sample and Abell clusters 
cross-correlation. Only accurate GRBs weighted by their positional errors were used for the 
evaluation. The shaded area shows the 5% probability for random GRBs sample to produce 
higher cross-correlation signal than the true GRBs sample signal. 

correlation alone, we get a rough estimate for the error multiplicative factor (e.g. Olivier et 
al. 1990) by letting E be the cluster mean surface density and w c i_ c i(6) ~ AO -1 . We then 
evaluate 



for R > 1 and here R > 0) and a single 4° radius bin, we obtain \J\ + u ~ 2.0. Though this 
estimate is only approximate, it allows a model-free interpretation of our results later on. 

In order to circumvent the difficulties in the cross-correlation evaluation, we have instead 
tested the null hypothesis, namely that the GRBs show no correlation with the rich clusters. 
We calculate N gr b_ c i(6) the number of GRBs— cluster pairs whose separation is smaller than 
a given angle 6 for 6 = 1°, 2°,..., 6°. We then create 500 random Poissonian realizations 
(artificial catalogs = "ac") of particles with the same selection function and error distribution 
as the GRBs and calculate N t ac _ d (6) for each realization i. We then attempt to rule out the 
null hypothesis by considering the fraction of random artificial catalogs for which N gr b_ c i(6) > 
N l ac _ d (9). This provides us with a direct estimate of the null hypothesis rejection level. 





to get Aw 



(1 + u) 1 / 2 (Aw)p. For expressed in degrees, A 0.7 (cf. Batuski et al. 1989 
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Figure 2: Rejection levels for the null hypothesis of no cross-correlation between GRBs and 
rich clusters. Different lines are for all GRBs and a sub-class of accurate position GRBs 
weighted or not by their positional errors. 

Since we conduct a purely statistical experiment we try to maximize the statistical 
significance by performing the same analysis for the full GRBs sample (\b\ > 30°) and for 
the "accurate" sub-class. In each of the above cases we have also carried out calculation 
for weighted correlation weighting each GRB - cluster pair by the GRB error, S~ 2 . Figure 
I shows the entire probability distribution for cross-correlation between artificial catalogs of 
the GRBs accurate weighted sample and the Abell clusters catalog. Only in 5% out of the 
500 artificial catalogs a correlation signal higher than the true correlation was found. The 
results for the different sub-samples and weighting schemes are summarized in table I and 
plotted in figure 2. 

Checking for self consistency within the table entries, we notice that for the last case 
("accurate" sub-class with statistical weighted cross-correlation) where we find the largest 
statistical significance, the maximum is obtained for 4° top-hat on the sky. The same holds 
for the second case ("accurate" sub-class with a regular correlation function) and for the 
third case (all bursts with weighted correlation function). These results reflect the combined 
effect of cross-correlation accompanied by systematic and random errors. Even if all GRBs 
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Table 1: Significance of cross-correlation coefficient (%) 
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reside in Abell clusters we could not expect higher correlation at smaller angles because of 
the errors. The statistical significance increases when we use weighted correlation function 
(column 5 vs. column 3) and when we use the accurate sub-class in place of the whole 
sample (column 5 vs. column 4 and column 3 vs. column 2). Both trends are reassuring. 
We recall that the more accurate bursts are stronger on average, hence nearer and therefore 
more correlated with the clusters. 

4 The Inferred Distances 

Having found a cross-correlation signal and the optimal angle to look for it, we can esti- 
mate the GRBs fraction that actually contribute to the signal, i.e. overlaps with the cluster 
sample volume. In a model where all GRBs reside in clusters, we compare the ratio between 
the detected cross-correlation and the two-dimensional top-hat average over the cluster auto- 
correlation. For a circle of radius 6 C we get (w c i_ c i(6 c )) ~ 2A6~ 1 and (tu c ;_ c ;(4 )} ~ 0.35. 
The measured, minimally biased values for the cross-correlation at a 4° bin are 0.013 and 
0.093 for all GRBs and for the "accurate" sub-sample respectively. While the former is 
not much of a use due to signal-to-noise ratio of ~ 1, the error estimate for the latter is 
Atu(4°) ~ 0.053. Notice that the statistical significance of the measurement has already 
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been established, based on the null hypothesis rejection, and the current analysis is needed 
for the overlap estimate only. The overlapping fraction is hence ~ 0.093/0.35 = 0.26 ± 0.15. 
For the "accurate" sub-sample with 136 GRBs, we conclude that 35 ± 20 of them lie in a 
volume of a 600 h _1 Mpc sphere. Using the relative amplitudes, and the expression for the 
expected angular cross-correlation (Lahav, Nemiroff, & Piran 1990, equation 9) one gets a 
characteristic depth for the GRBs "accurate" sub-sample of i?* j3r 6_ a = 940ti2o^ -1 Mpc i.e. 
Zg r b-a = 0.3lto!o4- I n a different estimates we calculate n(z), the fraction of observed bursts 
up to a redshift z (e.g. Cohen & Piran 1995) in a given cosmological model. Given the 
fact there are 35 ± 20 GRBs within z = 0.2, we obtain for = 1 and no source evolution 
z-max-a = 0.36io!o7 m which BATSE will detect 156 bursts. The equivalent estimate for all 
GRBs yields ~ 4% overlap with the cluster sample, namely ~ 22 GRBs within the same 
volume but with a large error. The much larger error can be interpreted as a steep build-up 
of the noise level due to non-overlapping GRBs. For the whole sample this implies a maxi- 
mum detection redshift of z max = 0.7to!o7- The two estimates are consistent with each other 
and suggest that the main correlation contribution for the entire GRBs sample essentially 
comes from the "accurate" sub-sample. An equivalent calculation, using the cluster-galaxy 
correlation function (Mo, Peacock, & Xia 1993; Lilje, & Efstathiou 1988), with approxi- 
mately half the cluster auto-correlation function amplitude in the relevant range and under 
the assumption that all GRBs reside in galaxies, leads to a larger overlapping fraction for 
the "accurate" sub-sample of ~ 0.5 and to lower z values. 

5 Conclusions and Discussion 

We have found a cross-correlation of GRBs with Abell clusters. The large rich cluster 
sample helps to overcome the statistical noise that prevents us from finding any GRBs auto- 
correlation. The correlation we found is believed to be a true spatial correlation. However, 
weak lensing magnification bias could also contribute to the observed effect. The predicted 
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magnification correlation averaged on a 6 C radius window is given by (e.g. Bartelmann 1995), 



where w^s is the matter-magnification correlation, b the considered lenses biasing factor and 
a is the power law index for the intrinsic GRBs source counts as a function of their flux. 
For the most favorable case for magnification we assume all GRBs are at z = 1.5, consider 
z < 0.2 clusters along the line of sight with a biasing factor 6 = 5, and use the standard 
CDM model function for these parameters (Bartelmann 1995) of w^s — 7 X 10~ 4 #~°' 6 . The 
obtained value is then (tu(4 )) < (—a — 1)0.014, which is one sixth of the detected signal for 
the accurate sample, and comparable to the (uncertain) signal for all GRBs (for reasonable 
values of a). 

If real spatial correlation is responsible for the detected signal, we know that Abell 
clusters extend up to ^ 0.15 — 0.2 and therefore the relative amplitudes of the cluster 
auto-correlation and the GRBs cluster cross-correlation provides an estimate for the GRBs 
redshift distribution. For a sub-sample of accurately measured GRBs we found that z max _ a 
of this population is 0.361q;q7 if all accurate GRBs are in clusters and z max _ a ~ 0.271q;q5 if 
they are all in galaxies. For the whole GRBs sample this implies z max = 0.7t.o'^ if GRBs 
originate in clusters and 0.5to!i5 if they are all in galaxies. The upper limits of these values 
are within the range of z max estimates as obtained from fitting the peak-flux distribution 
to a cosmological distribution. Cohen and Piran (1995) find z max = 2.1^;° for the long 
bursts and z max < 0.5 for the short bursts, while Fenimore et al. (1993) find 0.9 
for the combined population. However, In a more recent analysis Fenimore & Bloom (1995) 
find that if GRB time dilation is entirely attributed to cosmological redshift streching, then 
Zmax > 6. Here we have a combined population of long and short bursts. The average values 
are on the low side and this may indicate that lensing magnification and not just spatial 
correlation contributes to the observed signal. More data is needed to decrease the error 
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estimates and to determine which of the two possibilities is the dominant contribution. 

The cross-correlation is the first association of GRBs with any other population of 
astronomical objects. It demonstrates, of course, the cosmological origin of the bursts. It 
does not mean, however, that the GRBs originate necessarily in clusters, as it is consistent 
with the picture in which the bursts originate in galaxies or galaxy halos that in turn, are 
strongly correlated with Abell clusters. 

We thank Ehud Cohen, Ramesh Narayan, Reem Sari, and Eli Waxman for helpful 
discussions. This research was supported in part by the US National Science Foundation 
(PHY 91-06678) and by the Israeli National Science Foundation. 
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